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Correlated X-ray and optical variability in KV UMa 
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Abstract. We present coordinated X-ray and optical observations of the X-ray transient KV UMa (= XTE 
J1118+480) at high time resolution. The optical variation associated with spikes in the X-rays consist of a dip 
preceding the spike by 0-5 s, and a peak following it. The cross-correlation function between X-rays and optical has 
the same shape. It rises sharply, within 30ms of the peak of the X-rays. The shape of the cross correlation shows 
significant intrinsic variability on time scales as short as 30s. Analyzing this variability with principal component 
analysis (calibrated with simulated data), we find two statistically independent components. The first is similar 
to the average cross-correlation. The second also consists of a dip followed by a peak, but on a 3x shorter time 
scale. These properties of the optical light, together with the high optical brightness of the transient, are not 
easily explained by reprocessing of the X-rays. 
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1. Introduction 

The transient X-ray source XTE J1118+480 (=KV UMa) 
appeared in January 2000 (Remillard et al. 2000) and 
(with an interruption in February) was active as an X- 
ray source till end of July. The source appears to be a 
relatively nearby (2 kpc) binary with an orbital period 
of 4 hrs, and a K-star companion. Its mass function of 6 
M Q makes the primary a solid black hole candidate, (Mc 
Clintock et al. 2001a, Wagner et al. 2001). The X-ray out- 
burst was somewhat unusual compared with 'standard' 
transients such as X-ray Nova Muscae. The X-ray flux 
was relatively constant during the main outburst episode 
(March- July). The X-ray luminosity (1-100 keV) was only 
about 10 36 erg/s, and the source was always found in a 
X-ray hard state (most of the X-ray energy flux appearing 
around lOOkeV). Since hard state spectra are the rule dur- 
ing the decay phase of X-ray transients, this suggests that 
the accretion rate of the system was significantly below 
the Eddington value. 

Due to the low galactic absorption at the position of 
the source, its spectral energy distribution (SED) could be 
determined from the IR to X-rays, including the crucial 
EUV part of the spectrum (Hynes et al. 2000, McClintock 
et al. 2000b), a unique result for X-ray transients. These 
observations showed that the source had an unusually high 
optical-to- X-ray flux ratio compared with other transients. 
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The SED appears to peak at 100 eV, and the inferred UV- 
to-EUV flux of 10 36 erg/s was about the same as the X-ray 
flux. The optical-to-X-ray flux ratio was about 0.1. High 
optical brightness during a hard X-ray state has been seen 
before in the unusual long-term variable source GX 339-4 
(Motch et al. 1983). 

The long duration of this transient, its location above 
the galactic plane, and its relatively close distance made it 
an ideal target for simultaneous X-ray and optical obser- 
vations. Though such coordinated observations have been 
attempted regularly, successful observations at high (sub- 
second) time resolution have been reported previously in 
only two cases. Motch et al. (1981) reported correlated 
fast optical and X-ray variability from 96s of observations 
of GX 339-4, when it was in one of its rapidly variable and 
optically bright states. The X-ray-optical cross-correlation 
of these data was reported to show a negative correlation 
preceding the X-rays by a few seconds. Hynes et al. (2000) 
reported simultaneous HST and XTE observations of the 
black hole transient GRO J1655-40. The cross-correlation 
showed, in contrast with the case of 339-4, a positive cor- 
relation peak following the X-rays, interpreted by the au- 
thors as evidence of reprocessing of X-rays by the accretion 
disk. 

The optical emission of black hole candidates in out- 
burst thus appears puzzling, both in its unpredictable am- 
plitude, and its contradictory correlations with the X-rays. 
Theoretically, the optical light is interesting because it 
may provide a way to probe the structure of the accretion 
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flow, if the light is caused by reprocessing of X-rays from 
the central object (e.g. Arons & King 1993). Evidence for 
such reprocessing of X-rays by an accretion disk has been 
found in type-1 X-ray bursts from accreting neutron stars 
(Matsuoka et al. 1984, Turner et al. 1985, van Paradijs et 
al. 1990, Kong et al. 2000). 

On the other hand, lower limits to the optical bright- 
ness temperatures inferred from the short optical variabil- 
ity in GX 339-4 led Fabian et al. (1982) to conclude that 
a more likely mechanism for the optical radiation in the 
Motch et al. data is thermal cyclosynchrotron emission in 
a strong magnetic field. We reached a similar conclusion 
in a first description and interpretation of the results re- 
ported here (Kanbach et al. 2001). If this interpretation is 
correct, the regions contributing most to the optical-UV 
light would be much closer to the black hole than in a 
reprocessing model, at least in objects like KV UMa and 
GX 339-4, and optical light of black hole transients would 
potentially be a powerful new diagnostic of the inner ac- 
cretion flow. 

In Kanbach et al. (2001) we found that the average 
X-ray/optical cross correlation function consists of a peak 
starting at zero lag, preceded by a 'precognition dip'. It 
thus combines both of the effects seen previously by Motch 
et al. and Hynes et al. This raises the question how these 
two components of the cross correlation are related, in 
particular since the dip at negative lag is not very easily 
explained in models proposed so far. In this paper we ad- 
dress this question by a more detailed statistical analysis 
of our observations of KV UMa. The theoretical interpre- 
tation will be developed further in a future paper. 



2. Observations 

X-ray and optical observations were collected on the nights 
of 4,5,6 and 7 July 2001. X-ray observations were made by 
the Rossi X-ray Timing Explorer, for about one orbit on 
each of these nights. The total duration of useful simulta- 
neous observations was 2.5 hours. The optical observations 
were done with the OPTIMA photometer (Straubmeier et 
al. 200 1) attached to the 1.3 m telescope on Mt. Skinakas, 
Crete ( |http: / /observatory. physics. uoc.gr/| ), which is op- 
erated jointly by the Foundation for Research and 
Technology Hellas (FORTH) and the Max-Planck-Institut 
fur Extraterrestrische Physik. The photometer consists of 
a cluster of fiber-coupled avalanche diodes. The diodes are 
sensitive from 450 to 950 nm at a mean quantum efficiency 
of 50%. Individual photon arrival times are recorded to an 
accuracy of 2/xs with a GPS-based clock, at a maximum 
count rate of 200 000/s. On several of the nights, the tele- 
scope was facing into a fairly strong wind, causing mod- 
ulations of the signal at frequencies around 9 and 40 Hz 
due to telescope vibration. Most of the results reported 
here are based on the correlation between optical and X- 
ray signals, for which this artefact has little consequence. 
In the cross-correlation functions it appears as a source of 
noise but has no systematic effects. 




Fig. 1. Two examples of the X-ray and optical light 
curves, around the time of peaks in the X-rays. 
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Fig. 2. Average X-ray and optical light curves of 100 
peak-aligned events like the ones shown in Fig. |]. 

The average total X-ray count rate was 500 s _1 (3 
detectors) , typical for the brightness of the source during 
most of its outburst. The optical brightness was around 
V=13, resulting in a count rate of about 24000 s _1 . While 
the X-ray flux was highly variable, with amplitudes typical 
of the hard states of black hole transients or Cyg X-l, the 
rms variability of the optical flux was only about 10% on 
time scales of 100ms-30s. 



3. Analysis 

An example of the X-ray and optical light curves is shown 
in Fig. [|. The data show the original count rates, without 
background or baseline subtractions. Peaks in the X-rays 
with durations of a few seconds as seen in this example 
were fairly typical. Visual inspection of such events sug- 
gests that the corresponding optical signal consists of a 
dip, preceding the X-rays by several seconds, followed by 
a positive response peaking a few seconds after the X-rays. 

In order to quantify this impression in a more objec- 
tive way, an automatic peak-finding algorithm was devised 
to collect and superimpose these events. The X-ray time 
series was compared with a box-car smoothed version of 
itself (with a width of 30s). Time intervals where the ac- 
tual signal was more than twice the smoothed signal were 
identified as peaks. The 100 events found in this way were 
aligned in time on the maxima of the X-ray peaks. 

The average X-ray and optical light curves of the peak 
events found in this way is shown in Fig. B. This aver- 
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Fig. 3. Average X-ray/optical cross correlation from 297 sample of the time series. Right panel shows the cross 
correlation on an expanded scale; note the sharp rise at zero lag. 



age confirms the subjective impression given by Fig. [y, 
The dip in the optical preceding the X-rays appears to be 
a systematic effect, though its amplitude varies between 
events. 

3.1. Cross-correlation 

The superposition of events shown in Fig. |^ is limited in 
its time resolution by the width of the X-ray peaks, of the 
order of a second. A view at much higher time resolution, 
at the expense of a somewhat less direct interpretation, is 
given by the cross correlation between the data. The cross- 
correlation of the X-ray time series >Sx and the optical 
series So, CxoW — Sx* So, where a * denotes the con- 
volution, was computed for 297 independent segments of 
data of 25s length. The average of these cross-correlations 
is shown in Fig. |^. 

The cross correlation contains a residual of the pho- 
ton noise in the original time series. In Figs. ||, [| this 
shows up as the weak high-frequency noise component. 
Its amplitude is much lower, however, than the noise due 
to the intrinsic variability of the cross correlation, as Fig. 
[| shows. 

The cross-correlation appears to contain significant 
signal for lags (optical later than X-rays) from about - 
5 to +5s. On these time scales, the correlation shows the 
same properties as the direct comparison of the signals 
shown in Fig. ||: a dip followed by a peak of about the 
same amplitude. 

As the right panel of Fig. ||| shows, the change from 
the dip to the peak takes place quite close to zero lag. 
The slope of the correlation increases by a factor 10 within 
about 30 ms at t = 0, and 50% of the peak is reached at 
t = 100 ms. 

3.2. Dip and peak 

A peak starting at t — is understandable in a reprocess- 
ing model for the optical light, and its duration is also of 
the order expected for reprocessing of X-rays by an ac- 
cretion disk. Evidence for such reprocessing has been seen 



in X-ray bursters (Matsuoka et al. 1984, van Paradijs et 
al. 1990, Kong et al. 2000). It has also been advanced to 
explain the optical response in some black hole transients 
(O'Brien & Home 2000, Hynes et al. 1998). A delayed op- 
tical response is also a natural result in the outflow model 
of Kanbach et al. (2001). In this model, the unusually large 
optical flux of KV UMa does not come from the accretion 
disk but instead is produced by synchrotron emission in a 
magnetic outflow. The dip in the optical light preceding 
the X-rays is not very natural in either of these explana- 
tions, however. This raises the question whether the dip 
should be regarded as an independent phenomenon, or as 
a part of the correlation function associated directly with 
the peak. 

The question whether the dip represents an indepen- 
dent physical process from that producing the peak can be 
answered to some extent by investigating the variability 
of the dip and the peak between different samples of the 
cross correlation. If their amplitudes turn out to be statis- 
tically uncorrelated, it is unlikely that they are due to the 
same process. Significant variability is in fact present in 
both the dip and the peak. This is shown in Fig. ||, which 
presents a random subset of 16 of the 300 correlations. 
From this figure it is also clear, however, that a large part 
of the variability of the cross correlation appears to be in- 
dependent of both the dip and the peak. In many cases, 
this additional component masks the dip and/or the peak. 
From simple inspection of figures like these, it is hard to 
determine if this additional variability is just some form 
of intrinsic random variability in the source, or whether it 
contains additional systematic variations. 

4. Principal components of the cross correlation 

A more quantitative way of investigating the variability of 
a signal like the cross correlations in Fig. ^ is by a prin- 
cipal component analysis (PCA, e.g. Kendall 1980, not 
to be confused here with the PCA detectors on board of 
XTE). If the signal is variable, and a sufficiently large set 
of samples of it available, this method can decompose the 
samples into components whose amplitudes vary statisti- 
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Fig. 4. 16 randomly selected samples from the set of 297 X-ray/optical cross correlation functions. 



cally independently through the set. The method has been 
used in astronomy, among other applications, for the clas- 
sification of galaxy spectra from large surveys (e.g. Folkcs 
et al. 1999), AGN spectra (Boroson 2002), Chepheid light 
curves (Kanbur et al. 2002), and velocity distributions in 
molecular clouds (Brunt & Kerton 2002). The components 
found in this way can often be interpreted more easily in 
terms of variations in physically quantities than the orig- 
inal (large) datasets. 



In practice, significant care has to be excercised in in- 
terpretating the principal components. In the following we 
do this by comparison with forward modeling. First we ap- 
ply PCA to our set of 297 samples of the cross correlation 
and interpret a few of the most significant components. 
We then verify the stability and significance of these com- 
ponents by constructing sets of synthetic cross correlation 
functions, and processing them in the same way. These 
synthetic samples consist of superpositions of a few in- 
put components whose shapes are guessed, with statistical 
properties similar to those of the real set of cross corre- 
lations. The input components and the statistics of their 
amplitudes are adjusted to find the minimal number of 
'true' components, and their shapes, that are needed to 



reproduce the pattern of components found by the PCA 
of the real data set. 



4.1. Method 

For our application, the PCA is used as follows. Let the 
X-0 cross correlation function of sample j be denoted 
by aP(i). Let m be the number of available samples of 
the correlation function (m = 297 in the present ap- 
plication), discretize the time-delay axis t in n bins ti 
(—12.5 < ti < 12.5s, n = 150 in our case), and write 
a\ = &>{ti). We treat the n time bins as the variables 
whose statistical properties are to be found by PCA. Since 
neighboring time bins arc highly correlated, the number 
of degrees of freedom is much less than n. The degree of 
correlation of the time bins with each other contains in- 
formation on the actual number of degrees of freedom in 
the signal. We measure the correlation between time bins 
i and k by the covariance matrix Cik- 



Cik = ^2(a 3 t - - Ofc), 



(1) 



where a = ^\ a? /m is the average of a in the set. The 
principal components of a are then conventionally defined 
as those linear combinations of the columns of a] that 
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diagonalize the covariance matrix. [In many PCA appli- 
cations, the correlation matrix is used instead of the co- 
variance matrix. This is important if the variables mea- 
sured are different physical quantities. For our case, this 
makes practically no difference, since the variables a» are 
the discretized values of a single continuous function.] If 
the coefficients of these linear combinations are denoted 
by the matrix Cjfe, the principal components pk are given 

by 



p1 vor=0.380 av/dev=1.20 



p2 var = 0.285 av/dev = 0.24 



(2) 



and the samples themselves are given, in terms of the prin- 
cipal components, by the inverse of c, 



k 



(3) 



These components are somewhat abstract, since a lin- 
ear combination of different time bins is not a physically 
very meaningful quantity in our application. A more use- 
ful representation of the principal components is found by 
noting that the quantity 



j —i j 

a ki — c ik Pk 



(4) 



describes the amplitude of principal component k in the 
i-th time bin of the j-th sample. The average of this am- 
plitude over the set of samples, 



3 



(5) 



describes the principal components in the 'time domain', 
i.e. as a decomposition of our average cross correlation 
function. We use this representation for the rest of our 
analysis. The average of the cross correlation function over 
the set, a, is the sum of these principal components: 



Eft 



(G) 



4.2. Results 

The covariance matrix is singular if the number of time 
bins n is larger than the number of samples m. This limits 
the time resolution that can be used for the discretization 
of our cross correlation functions to 300 points^. The re- 
sults depend somewhat on the range in t used for the 
analysis. Most of the signal of interest appears between 
t = ±5s. Using this range, and bins of 30ms yields the 
principal components shown in Fig. ^|. 

1 This limitation can be circumvented by data-compression 
of the samples. A convenient compression is obtained by apply- 
ing a wavelet-transform to the samples. PCA is then applied 
on the first no < m wavelet coefficients, and the resulting prin- 
cipal components transformed back to the time domain by the 
inverse wavelet transform. For no — n this process is loss-free 
and yields identical results 
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p5 var = 0.034 av/dev = 0.45 



p6 var-0.026 av/dev = 0.01 
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Fig. 5. The first six principal components from a PCA of 
the set of 297 samples of the X-ray/optical cross correla- 
tion. Numbers above the panels show the fraction of the 
variance in the cross correlation accounted for by the com- 
ponents, followed by ratio of the average to the deviation 
of the component's amplitude. 



The amplitude A of components varies both in abso- 
lute value and sign. For the first component we find that 
the average is larger than the variation, but for the higher 
components the average is generally small compared to 
the variation. Instead of showing the components at their 
average amplitude A, we have therefore plotted them at 
a scale that takes into account both the average and the 
variation. The value used for this, somewhat arbitrarily, 
is A = {A 2 + a\) x / 2 , where a a is the rms deviation of A. 

The first two components together account for 67% 
of the variance between the samples of the cross corre- 
lation function, if a time-lag interval (—At < t < At) 
with At = 5 s is used. The first component appears quite 
similar to the average cross correlation. The second com- 
ponent peaks around the 'precognition dip', though it has 
the opposite sign and extends somewhat to positive lag 
as well. If the time interval used is shorter, the first two 
components account for a larger fraction of the variation, 
and less for a longer interval. For At = 2s the first two 
components account for 86% of the variance, for At = 12 
s only 44%. This is understandable if the cross correlation 
consists of a signal of relatively stable shape concentrated 
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Fig. 6. Stability of the pincipal components with respect 
to time variability of the data. The full set of correla- 
tion functions has been divided into halves in two ways. 
The dark lines show the even and odd numbered samples 
(chronologically ordered), measuring differences on a time 
scale of 25 s. The lighter lines show the first and second 
halves of the set, indicating differences on a time scale of 
hours to days. The higher components appear to be quite 
stable in this sense, but in fact most of these turn out to 
be artefacts (see text and Figs. |j). 



around t = 0, plus an unsystematically varying compo- 
nent that is spread out over a wider range in the time 
delay coordinate t. 

The nature of the remaining components is not obvious 
at this point. Some of them, such as P3 and P5, appear 
to contain signals related to the sharp rise of the cross 
correlation at t = 0. The PCA also interprets pure noise 
in the data in terms of additional components, and it is 
not a priori obvious which components are systematic and 
which may be just artefacts. A useful test to see which 
of the components is stable with respect to noise is to 
divide the set of samples in two halves and comparing 
the components derived from each of these subsets. This 
is shown in Fig. ^|. The four curves shown result from 
two different ways of dividing the set in halves. In the 
first, the two subsets consist of the even numbered and the 
odd numbered samples of the (chronologically ordered) 
set. The difference between these two subsets measures the 
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Fig. 7. Top: input components for the synthetic set of 
cross correlations. Their amplitudes are varied indepen- 
dently. In addition there is a noise component (not shown) 
whose contribution to the variance peaks around 0.2 Hz. 
The average cross correlation produced by these compo- 
nents is shown in the bottom panel. 



intrinsic variability of the cross correlation on a time scale 
of 25s (the duration of a sample). In the second division, 
the first half of the set is compared with the second half. 
This measures the variability on time scales of days. 

As Fig. ^ shows, the components, especially the higher 
ones, are remarkably stable on both the short and the long 
time scales. From this one would be tempted to conclude 
that all six components shown are real. The process by 
which the components are derived in PCA, however, leaves 
the possibility open that some of the components are arte- 
facts. This is because the components are constructed as 
orthogonal directions in the abstract n— dimensional space 
of variables. The directions of insignificant higher com- 
ponents (constructed out of noise) are thus biased in a 
systematic way by the directions of the significant lower 
components. 

The only way to test which of the components is real in 
this sense is by comparison of the results with a synthetic 
data set. 

4.2.1. Tests with synthetic data 

To test the significance of the components, we have con- 
structed synthetic data and applied PCA on this data in 
the same way as was done with the real data. By vary- 
ing the input data, the resulting PCs can be adjusted for 
agreement. Since many of the components found are sus- 
pected to be spurious, the aim of the excercise was to find 
the minimum number of independently varying signals 
needed to reproduce the components found by PCA of the 
real data. It is not neccessary to start the synthetic pro- 
cess with time series of the count rates, since the questions 
to be answered concern the behavior of principal compo- 
nent analysis as applied to our cross-correlation functions. 
Thus, the synthetic data constructed consist of samples of 
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the cross correlation function. Each of the samples consists 
of a linear combination of a number of adjustable functions 
of the time delay coordinate t, here called 'signals'. The 
shape of these signals is the same in all samples, but their 
amplitudes differ. In addition to these signals, we add a 
random noise component that is different in each sample. 

The amplitudes of the signals are assigned randomly 
between the samples, normally distributed with specified 
means and deviations. The noise component consists of 
white noise, low-pass filtered by boxcar averaging with a 
width St. The parameters of the synthetic data are thus: 
the shapes of the signals, the means A\ and A^ and devi- 
ations <j\, (72 of their amplitudes, the amplitude A n of the 
noise component, and St. 
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Fig. 8. Components found by PCA of a sample of 297 
synthetic cross correlations, constructed from the two in- 
put components shown in Fig]?]. The shapes of the input 
components and their amplitudes were adjusted for agree- 
ment of the resulting synthetic PCA components with the 
first 6 components of the actual data (Fig. ||) 



PCA on a set of 297 of these synthetic cross corre- 
lations yields the components shown in Fig. ||, for a case 
with two input signals plus a noise component. The shapes 
of the two input signals were adjusted to fit the shapes of 
the output components to those of the real data (Fig. 5). 
The final result for these input signals is shown in fig. 1 

The amplitude of first two PCs of the simulated data 
are related directly to the average amplitudes of the two 
input signals. The amplitudes of PCs 3 and higher, how- 
ever, turn out to be proportional to the amplitude of 
the boxcar-filtered noise component instead, showing that 
they are artefacts of the PCA. PCA somehow constructs 
these higher components out of the noise component. The 
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Fig. 9. Average power spectra of the X-ray and optical 
time series (unit on y-axis: variance per factor e in fre- 
quency). The peaks in the optical spectrum around 10 Hz 
are artefacts due to telescope vibrations. 



shapes of these higher components, however, depend sen- 
sitively on the shape of the two input signals. The input 
signals can be determined rather accurately by the fitting 
process. Two signals plus a noise component turned out 
to be the minimum number needed. Tests with only one 
signal plus noise have also been done. These are able to 
reproduce the first two PCA components, but leave most 
of the detail in the higher components unexplained. 

Except for the first, none of the components found by 
the PCA resembles the input signals very well. It is thus 
clear that PCA without comparison with a simulation us- 
ing synthetic data can not be used to identify the sta- 
tistically independent components reliably in the present 
problem. Unfortunately, the forward modeling process was 
quite cumbersome, owing to the large number of parame- 
ters to be adjusted, and is not easily automated. 

From this test we conclude that just two independent 
signal components are needed to explain the observed vari- 
ability of the cross correlation, apart from an aditional un- 
systematic noise component. The interpretation of these 
components presents some puzzles. One might have hoped 
that the actual independent signals would be the 'spike' 
and the 'precognition dip', but this turned out not to be 
the case. Extensive attempts were made to reproduce the 
PCA components with input signals resembling a pure dip 
before t = and a spike after t = 0, but no satisfactory 
fit was found. When the shape of the signals was relaxed 
away from this prejudice, iteration converged onto the 
signals shown in Fig. [7| Their shapes could be determined 
rather accurately by this fitting process. 

Both signal components consist of a dip followed by 
a peak. The second signal component resembles a com- 
pressed (in time) version of the first, but with the peak 
being somewhat less pronounced relative to the dip. A 
possible interpretation is thus that the intrinsic shape of 
the optical response is a combination of a dip followed by 
a sharp peak, while the time scale of this pattern varies, 
by a factor of 3 or so. 

5. Power spectra 

Information complementary to the cross correlations is 
given by the power spectra. The average X-ray and optical 
power spectra are shown in Fig. ^. The broad peak near 
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Fig. 10. Principal components of the variability in the X- 
ray power spectra. (Note that the vertical scale is linear, 
instead of the log scale in Fig. 0) 



0.1 Hz in the X-ray power spectrum is associated with a 
QPO around this frequency. This QPO has been found in 
most of the X-ray data from KV UMa (Revnivtsev et al. 
2000, Wood et al. 2000), and with high amplitude also in 
the optical light curve (Haswell et al. 2000, Pavlenko et al. 
2001). In our data its amplitude was significantly smaller 
than in these earlier observations, and undetectable in the 
optical. Its frequency varied a bit during our observations, 
with some preference for frequencies of 0.08 and 0.2 Hz. 

The variation of the power spectrum can be studied 
with dynamical power spectra; such studies have been re- 
ported for many sources (cf. van der Klis 2000). In the 
spirit of our analysis of the cross correlation, the vari- 
ability in the power spectra can also by studied by PCA. 
Treating the 297 individual power spectra in the same way 
as the cross correlations yields components shown in Fig. 

The first component of the X-ray power spectrum is 
rather similar to the average, with a broad peak around 
0.1 Hz. The second component shows a dip at 0.1 Hz and 
peaks at 0.08 Hz and 0.2 Hz. This could be explained by 
the presence of QPOs at frequencies of 0.08, 0.1 and 0.2 
Hz, with the amplitude of the 0.1 Hz peak decreasing when 
the others are increasing. Tests with synthetic data show, 
however, that PCA also produces such a pattern if there 
is just a single a peak at a somewhat variable frequency. 
This seems a more likely interpretation in view of what is 
known about QPOs in similar sources. 

The optical power spectra are affected by telescope 
vibrations at frequencies near 10 and 40 Hz. As explained 
above, these vibrations were caused by the rather strong 
wind loading of the telescope on several of the nights. The 
optical variability is also contaminated with seeing effects, 
which have frequencies in the range range of interest here. 
Since the X-0 cross correlation reaches maximum values 
around 40% it is unlikely that seeing effects dominate the 
power spectrum, but its precise shape can not be reliably 
deduced from our observations. 

With this caveat, the optical power spectrum shown 
in Fig. |^ is much broader than the X-ray power spectrum, 
and there is no clear evidence for the QPO seen in the 
X-rays. It thus seems unlikely that a QPO contributes 
strongly to the X-0 cross correlations. The effect of the 
X-ray QPO's on our cross correlations may well be just a 
a random noise contribution. 



6. Conclusions and discussion 

We have studied the correlation of X-ray and optical prop- 
erties of the black hole transient KV UMa. Thanks to the 
length of the dataset, the statistical properties of this cor- 
relation could be studied in detail. By analysis of the time 
series themselves, we find that the optical response to an 
X-ray peak consists of a broad dip followed by a sharper 
peak. By cross-correlating the two time series the same 
pattern is found, but with a much higher signal-to-noise 
ratio. From the cross correlation we find that the onset of 
the peaks in the optical response coincides with the X-ray 
maximum to within 30 ms, and their initial rise is very 
fast. The minimum of the dip in the optical precedes the 
X-rays by about 2 seconds. 

The agreement of the pattern of the cross-correlation 
with that found directly in the time series establishes, on 
the one hand, that the dip and the delayed peak are both 
properties of the optical light curve rather than the X-rays, 
and on the other hand that they are not due to artefacts of 
background subtraction or the cross-correlation process. 

The shape of the cross correlation was found to be 
highly variable on time scales from minutes to days. A 
principal component analysis (PCA) of this variablity was 
found to be of some value in studying its statistical prop- 
erties. It turned out that the principal components them- 
selves, except for the first, are not very useful since their 
shape is highly affected by artefacts arising from a strong, 
randomly variable additive noise component in the cross 
correlation. 

We find that the PCA components can be used, how- 
ever, as an intermediate representation of the data, which 
can then be interpreted by a forward modeling process. 
The result of such a simulation shows that there are two 
systematic components which vary statistically indepen- 
dently. Both have the shape of a broad dip followed by 
a sharper peak, but their duration differs by a factor of 
about 3. 

The most puzzling aspect of the optical variations is 
of course the 'precognition dips'. If the optical light is due 
to reprocessing of X-rays, then an optical signal preced- 
ing the X-rays would have to be due to something in the 
accretion flow preceding an X-ray event. In principle, this 
could happen if the increase in mass flux which somewhat 
later is observed as an increase of the X-ray flux, starts at 
some distance from the hole, such that the accretion time 
scale from that distance is of the order of a few seconds. 
For typical accretion disk models, this would be a distance 
intermediate between the X-ray emitting region and the 
optical reprocessing regions in the outer disk. If this event 
is associated with a local thickening of the accretion disk, 
and the X-ray emitting region is small enough, such a 
'bump' on the disk surface might be sufficient to shield 
the outer disk temporarily from the X-rays, leading to a 
dip in the optical emission. 

While this possibility may be attractive in the repro- 
cessing model, reprocessing itself has difficulties explain- 
ing the shape of the optical response. The sudden change 
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in slope at zero lag turns out to be impossible to reproduce 
with any positive definite disk reprocessing function. The 
main reason for this is that the sudden change in slope 
requires much more variability on very short time scales 
(< 100 ms) than the X-ray light curve actually has. A 
second difficulty is that the optical-to-X-ray flux ratio in 
KV UMa was quite high, or the order 10%. If, as the spec- 
tral energy distribution of McClintock et al. suggests, the 
UV-EUV emission is part of the same process producing 
the visible light, the reprocessed flux would have to be as 
large as the X-rays themselves. 

From these difficulties it appears that alternatives may 
have to be explored to explain the visible light in black 
hole accreters like KV UMa and GX 339-4. In a future 
paper we will explore the possibility that the visible light 
is produced as thermal cyclo-synchrotron emission in the 
inner 20000 km of the accretion disk, as first suggested by 
Fabian et al. (1982) for GX 339-4 (see also Merloni et al. 
2001 for an interpretation of the spectrum of KV UMa). If 
this interpretation is correct, it suggests a new and poten- 
tially very powerful means of probing the inner accretion 
flow by means of optical, and especiallly coordinated X- 
ray and optical observations. 
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